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ABSTRACT 

We extend the formalism for the calculation of the relativistic corrections to the Sunyaev- 
Zel'dovich effect for clusters of galaxies and include the multiple scattering effects. We present a 
systematic method for the inclusion of the multiple scattering effects. The multiple scattering 
contribution is found to be very small compared with the single scattering contribution. For 
high-temperature galaxy clusters of ksT e w 15kcV, the ratio of the both contributions is —0.3% 
in the Wien region. In the Raylcigh- Jeans region the ratio is —0.03%. Therefore the multiple 
scattering contribution is safely neglected for the observed galaxy clusters. 

Subject headings: cosmic microwave background — cosmology: theory — galaxies: clusters: 
general — radiation mechanisms: thermal — relativity 
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1. INTRODUCTION 

Compton scattering of the cosmic microwave background (CMB) radiation by hot intracluster gas 
- the Sunyaev-Zel'dovich effect (Zel'dovich & Sunyaev 1969; Sunyaev & Zel'dovich 1972, 1980a, 1980b, 
1981) — provides a useful method to measure the Hubble constant H (Gunn 1978; Silk & White 1978; 
Birkinshaw 1979; Cavaliere, Danese, & De Zotti 1979; Birkinshaw, Hughes, & Arnaud 1991; Birkinshaw 
& Hughes 1994; Myers et al. 1995; Herbig et al. 1995; Jones 1995; Markevitch et al. 1996; Holzapfel et 
al. 1997; Furuzawa et al. 1998; Komatsu et al. 1999). The original Sunyaev-Zel'dovich formula has been 
derived from a kinetic equation for the photon distribution function taking into account the Compton 
scattering by electrons: the Kompaneets equation (Kompaneets 1957; Weymann 1965). The original 
Kompaneets equation has been derived with a nonrelativistic approximation for the electron. However, 
recent X-ray observations have revealed the existence of many high-temperature galaxy clusters (David 
et al. 1993; Arnaud et al. 1994; Markevitch et al. 1994; Markevitch et al. 1996; Holzapfel et al. 1997; 
Mushotzky & Scharf 1997; Markevitch 1998). In particular, Tucker et al. (1998) reported the discovery 
of a galaxy cluster with the electron temperature ksT e = 17 .4 ± 2.5 keV. Rephaeli and his collaborator 
(Rephaeli 1995; Rephaeli & Yankovitch 1997) have emphasized the need to take into account the relativistic 
corrections to the Sunyaev-Zel'dovich effect for clusters of galaxies. 

In recent years remarkable progress has been achieved in the theoretical studies of the relativistic 
corrections to the Sunyaev-Zel'dovich effects for clusters of galaxies. Stebbins (1997) generalized the 
Kompaneets equation. Itoh, Kohyama, & Nozawa (1998) have adopted a relativistically covariant formalism 
to describe the Compton scattering process (Berestetskii, Lifshitz, & Pitacvskii 1982; Buchlcr & Yuch 
1976), thereby obtaining higher-order relativistic corrections to the thermal Sunyaev-Zel'dovich effect in 
the form of the Fokkcr-Planck expansion. In their derivation, the scheme to conserve the photon number at 
every stage of the expansion which has been proposed by Challinor & Lasenby (1998) played an essential 
role. The results of Challinor & Lasenby (1998) are in agreement with those of Itoh, Kohyama, & Nozawa 
(1998). The latter results include higher-order expansions. Itoh, Kohyama, & Nozawa (1998) have also 
calculated the collision integral of the Boltzmann equation numerically and have compared the results with 
those obtained by the Fokker-Planck expansion method. They have confirmed that the Fokker-Planck 
expansion method gives an excellent result for ksT e < 15keV, where T e is the electron temperature. For 
ksTe > 15kcV, however, the Fokker-Planck expansion results show nonnegligible deviations from the results 
obtained by the numerical integration of the collision term of the Boltzmann equation. 
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Nozawa, Itoh, & Kohyama (1998b) have extended their method to the case where the galaxy cluster 
is moving with a peculiar velocity with respect to CMB. They have thereby obtained the relativistic 
corrections to the kinematical Sunyaev-Zel'dovich effect. Challinor & Lasenby (1999) have confirmed the 
correctness of the result obtained by Nozawa, Itoh, & Kohyama (1998b). Sazonov & Sunyaev (1998a, 
b) have calculated the kinematical Sunyaev-Zel'dovich effect by a different method. Their results are in 
agreement with those of Nozawa, Itoh, & Kohyama (1998b). The latter authors have given the results of 
the higher-order expansions. 

Itoh, Nozawa, & Kohyama (2000) have also applied their method to the calculation of the relativistic 
corrections to the polarization Sunyaev-Zel'dovich effect (Sunyaev & Zel'dovich 1980b, 1981). They 
have thereby confirmed the result of Challinor, Ford, & Lasenby (1999) which has been obtained with a 
completely different method. Recent works on the polarization Sunyaev-Zel'dovich effect include Audit & 
Simons (1999), Hansen & Lilje (1999), and Sazonov & Sunyaev (1999). 

As stated above, Itoh, Kohyama, & Nozawa (1998) have carried out the numerical integration of the 
collision term of the Boltzmann equation. This method produces the exact results without the power series 
expansion approximation. Sazonov & Sunyaev (1998a, b) have reported the results of the Monte Carlo 
calculations on the relativistic corrections to the Sunyaev-Zel'dovich effect. In Sazonov & Sunyaev (1998b), 
a numerical table which summarizes the results of the Monte Carlo calculations has been presented. This 
table is of great value when one wishes to calculate the relativistic corrections to the Sunyaev-Zel'dovich 
effect for galaxy clusters of extremely high temperatures. Accurate analytic fitting formulae would be still 
more convenient to use for the observers who wish to analyze the galaxy clusters with extremely high 
temperatures. For this purpose, Nozawa et al. (2000) have presented an accurate analytical fitting formula 
of 0.1% accuracy for the numerical results for the relativistic corrections to the thermal Sunyaev-Zel'dovich 
effect for clusters of galaxies. For the analyses of the galaxy clusters with extremely high temperatures, 
the results of the calculation of the relativistic thermal bremsstrahlung Gaunt factor (Nozawa, Itoh, & 
Kohyama 1998a) and their accurate analytic fitting formulae (Itoh et al. 2000) will be useful. 

In this series of papers devoted to the study of the relativisitc corrections to the Sunyaev-Zel'dovich 
effect for clusters of galaxies, we have so far restricted ourselves to the case of single Compton scattering. 
This is justified because the optical depth for the Compton scattering of the CMB photon inside the 
galaxy clusters is generally about 10~ 2 or smaller (Birkinshaw 1999). Nevertheless, it would be desirable 
to evaluate the effects of the multiple Compton scattering of the CMB photon inside the galaxy clusters 
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accurately, as we have already developed the method to calculate the relativistic corrections to the 
Sunyaev-Zel'dovich effect for the galaxy clusters with high accuracy. The multiple scattering effects have 
been already considered by many authors (see Birkinshaw 1999 for references). In this paper we wish to 
evaluate the multiple scattering effects in the same theoretical framework of this series of papers. 

The present paper is organized as follows. In § 2 we give the method of the calculation and the results. 
In § 3 we give discussion of the results and concluding remarks. 



MULTIPLE SCATTERING CONTRIBUTION 



In the present paper, we would like to derive the analytic expression for the multiple scattering 
contribution in the Sunyaev-Zeldovich effect for clusters of galaxies. As a reference system, we choose the 
system that is fixed to the center of mass of the galaxy cluster. The galaxy cluster is assumed to be fixed to 
the cosmic microwave background (CMB). Following Itoh, Kohyama & Nozawa (1998), we start with the 
Fokker-Plank expansion for the time evolution equation of the CMB photon distribution function n(ui): 
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In equation (2.4), W is the transition probability of the Compton scattering, ,f(E) is the relativistic 
Maxwellian distribution function for electrons with temperature T e . We have integrated equation (2.4) 
analytically with power series expansions of the integrand. The expansion parameter is 

k B T e 
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The explicit forms for Ik are given in Itoh, Kohyama & Nozawa (1998). 

We first assume the initial photon distribution of the CMB radiation to be Planckian with temperature 



To: 



n(X) = n Q (X) = (2.6) 



where 



Assuming T a /T e <C 1, one obtains the following expression for the fractional distortion of the photon 
spectrum derived by Itoh, Kohyama & Nozawa (1998): 

A„(.Y) y^Xe*_ [Yo + 0eYl + e 2y 2 + 9 3y 3 + ^ , (2 . 8) 
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y = a T f dhNeih), (2.9) 
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where <jt is the Thomson scattering cross section, N e is the electron number density, and the integral 
is over the photon path length in the cluster. The explicit forms for Yo,Yi, Y2,Y^ and Y4 are given in 
Itoh, Kohyama & Nozawa (1998). Equation (2.8) is the single scattering contribution, i.e. the first-order 
term in y. If the cluster of galaxies is optically thin, i.e. y <C 1, the single scattering approximation is a 
good approximation. In fact, the approximation is valid for most of the clusters. However, it is extremely 
important to calculate the next-order contribution in order to obtain more accurate theoretical prediction 
for the future observation of the Sunyaev-Zeldovich effect for clusters of galaxies. 

We now calculate the multiple scattering contribution. Since y <C 1 is realized for most of clusters of 
galaxies, the second-order contribution is considered to be sufficient. We now assume that the initial photon 
distribution has a first-order perturbation. Namely, 

n(X) = m{X) = n Q (X) + An(X), 

An(X)\ 
n 



noPOU + T-^j. ( 2 - 10 ) 



where the second term in equation (2.10) is given by equation (2.8). Inserting equation (2.10) into RHS 
of equation (2.1), and performing the standard calculation, we obtain the following expression for the 
fractional distortion of the photon distribution function including the second-order contribution: 
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In equation (2.11), the first term corresponds to the first-odrder contribution and the second term 
corresponds to the second-order contribution of the multiple scattering. In deriving equation (2.11), we 
have used the following identity relation for y: 



= \y 2 - (2.17) 

We have also neglected terms higher than O{0 4 ) in the y 2 contributions in equation (2.11). It is important 
to note that equation (2.11) satisfies the photon number conservation. With equation (2.11), we define the 
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distortion of the spectral intensity as follows: 

A, . ^^ffi = A/,+A/,. (2 , 8) 
e x - 1 n (X) 

The first term Aii contains a factor y, whereas the second term AI 2 contains a factor y 2 . In Figure 1 we 
show AI 2 /y 2 as a function of X for the case fcsT e = lOkcV. It is clear that the magnitude of AI 2 /y 2 has a 
maximum value atl~5 for /csT e =10keV. 

In order to estimate the relative importance of the multiple scattering contribution, we now define the 
following ratio: 

r = (2.19) 

In Figure 2 we show T for X = 5 as a function of the electron temperature T e . It is clear that T increases 
with a negative sign as the temperature of the cluster of galaxies increases, i.e., T w —0.3 at ksT e = 15keV. 
However, the multiple scattering contribution is small because of a further factor y. Namely, for the cluster 
of galaxies of k B T e = 15keV, we have 

AIn 

— - = yT «-0.3y w-0.3%, (2.20) 

where we used a typical value y w 0.01 of the galaxy clusters. Therefore the maximum effect of the multiple 
scattering contribution is —0.3% of the single scattering contribution for the observed high-temperature 
galaxy clusters. 

In the Rayleigh- Jeans limit where X — > 0, equation (2.11) is further simplified: 

2v q (i 17 o + 123 fl 2 1989 3 14403 4 

+2y 2 e 2 e {i- l le e + ^e 2 ^ . (2.21) 

With equation (2.21) we have the multiple scattering contribution for ksT e =15keV and y = 0.01 as follows: 

AIo 

— - w -y0 e k -0.03%. (2.22) 
Aii 

In the Rayleigh- Jeans region the multiple scattering contribution is safely neglected. 



3. DISCUSSION AND CONCLUDING REMARKS 



From the results presented in the previous section it is clear that the multiple scattering contribution 
AI 2 is very small compared with the single scattering contribution AI\. For high-temperature galaxy 
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clusters of ksT e ~ 15keV, we obtain the ratio A/ 2 /A/i w —0.3% at X = 5. In the Rayleigh- Jeans region we 
have A/2/A/1 « —0.03%. Therefore it is concluded that the multiple scattering contribution to the thermal 
Sunyaev-Zel'dovich effect for galaxy clusters can be safely neglected. The reader is therefore referred to the 
previous four papers in this series of papers which deal with the single scattering contribution in detail. 

This work is financially supported in part by the Grant-in-Aid of Japanese Ministry of Education, 
Science, Sports, and Culture under the contract #10640289. 
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Figure Captions 



Figure 1. The spectral intensity distortion AI 2 /y 2 as a function of X for the case /csT e = lOkeV. The 
dotted curve shows the contribution up to Z - The dashed curve shows the contribution up to 9 e Z\. 
The solid curve shows the full contribution up to 6^Z 2 . 

Figure 2. The ratio T as a function of fcsT e for a fixed value of X = 5. The dotted curve shows the 
contribution up to Z . The dashed curve shows the contribution up to 6 e Z\. The solid curve shows 
the full contribution up to 6\Zi. 



